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Abstract. In this paper we analyze the behavior of the Laplace operator with Neumann boundary con- 
ditions in a thin domain of the type ij" = {{xi,X2) £ \ xi £ (0, 1), < X2 < eG{xi, xi/e)} where the 
function G{x,y) is periodic in y of period L. Observe that the upper boundary of the thin domain presents 
a highly oscillatory behavior and, moreover, the height of the thin domain, the amplitude and period of the 
oscillations are all of the same order, given by the small parameter e. 

Resume: Dans cet article, nous analysons le comportement de I'operateur de Laplace avec conditions aux 
limites de Neumann dans un domaine fine du type ij" = {(xi, X2) £ | xi £ (0, 1), < X2 < eG{xi,xi /e)} 
lorsque la fonction G{x, y) est periodique dans la variable y de periode L. On observe que la limite superieure 
du domaine fine presente une comportement hautement oscillatoire et, en outre, 1' hauteur du domaine, 
I'amplitude et la periode des oscillations sont tons du meme ordre, donne par un petit parametre e. 

Keywords: Thin domain, oscillatory boundary, homogenization. 



In this work we study the asymptotic behavior of the sohitions of the Neumann problem for the Laplace 
operator 



where Ge(-) is a function satisfying < Gq < Ge(-) < Gi for some fixed positive constants Go, Gi and such 
that oscillates as the parameter e — > 0. We may think, for instance, that the function Gc is of the form 
Gf{x) = a{x) + b{x)g{x/e") where a, 6 : / 1— M are piecewise C^-functions defined on / = (0, 1), g : IR — ?► M is 
an L-periodic smooth function and a > 0, see Figure 1. This includes the case where the function Ge(-) is a 
purely periodic function, for instance, G^{x) = 2 + sin(a;/e") but also includes the case where the function 
Ge is not periodic and the amplitude is modulated by a function. As a matter of fact, we will be able to 
treat more general cases, see hypothesis (H) below, but to stay the general ideas in the introduction we may 
consider the proptotype function Gt{x) = a{x) + b{x)g{x/e°'). 

The existence and uniqueness of solutions for problem for each e > is garanteed by Lax-Milgram 
Theorem. We are interested here in analyzing the behavior of the solutions as e — J' 0, that is, as the domain 
gets thinner and thinner although with a high oscillating boundary. 

Observe that the domain is thin since C (0, 1) x (0, eGi) and its upper boundary oscillates due to the 
function ^(x/e"), (if a > and g is not a constant function). 
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1. Introduction 




(1.1) 



with € L'^{R'^), j/*^ = (i/|, 1^2) is the unit outward normal to dR'^ and is the outside normal derivative. 
The domain R"^ is a thin domain with a highly oscillating boundary which is given by 

R" = {{xuX2) G I e (0,1), 0<X2< eG,{xi)} 
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Figure 1. The thin domain 



Since the domain is thin, "approaching" the hne segment (0, 1) C M, it is reasonable to expect that 
the family of solutions w'^ will converge to a function of just one variable and that this function will satisfy 



an equation of the same type of (1.1) but in one dimension, say Lu + u = /i in (0, 1) with some boundary 
conditions, where L is a second order elliptic operator in one dimension. As a matter of fact, if the function 
Ge(x) is independent of e, (say a = or g = 0), that is, the thin domain does not present any oscillations 
whatsoever: 

R" = {{xi,X2) : Q < xi <l,Q < X2 < e.G{x)} 

the limit equation is given by 

-^{G{x)w,),+w{x)^ f{x) xe(0,l) 
wM = = 

see for instance [101 E] . Observe that the geometry of the thin domain enters into the limit equation through 
the diffusion coefficient. 

Moreover, if we consider < a < 1 and if we assume that a{x) + b{x)g{x/e°') h{x) w-L^(0, 1) and 
a{x)+b(x)g(x/t'') ^^^) w-i^(0, 1) (obscrve that h{x)k{x) > 1 a.e. and in general it is not true that 
h{x)k{x) = 1), then the limit problem is 

see [2. Observe that this case contains the previous o ne. If a = 0, then h{x) ~ a(x) + b{x)g{x) = G{x) and 
Hx) = a{a:]+b{x)g{x) = ' ^'''^ rccover equation Ol. 



In this work we are interested in addressing the case a = 1, that is Ge{x) = a{x) +b{x)g{x/e), where none 
of the techniques used to solve the previous ones apply. Observe that this situation is very resonant: the 
height of the domain, the amplitude of the oscillations at the boundary and the period of the oscillations are 
of the same order e. Moreover we are interested in addressing not only the purely periodic case, that is, the 
case where the function G^{x) = G{x/e) for some L-periodic smooth function G but also the general case 
where the amplitude of the oscillation depend on a; in a continuous fashion, that is, in our prototype case, 
the situation where a and b are not piecewise constant, but piecewise continuous function. 

The purely periodic case can be addressed by somehow standard techniques in homogcnization theory, as 
developed in [SI [HI HI]. If G^{x) = G{x/e) where G is an L-periodic function and if we denote by 

y* = {(yi, 2/2) e : < t/1 < L, < 2/2 < G{vi)} 

then the limit equation is shown to be 

i -qoWxx+w = ,f{x), xG(0,1) 
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where 



90 = 



and X is the unique solution of 



-AX = in y* 
dX 

X{0,y2) = X{L,y2) on Bq 



X dyidy2 = 



(1.4) 



where Bq is the lateral boundary, Bi is the upper boundary and B2 is the lower boundary of dY* , that is 

Bo = {(0, y2) : < 2/2 < G(0)} U {(L, 2/2) : < zj2 < G{L)} 
= {(2;i,G(yi)) :0<yi <L} 
B2 = {iyi,0):0<yi<L}. 

We refer to for a complete analysis of the purely periodic case of a nonlinear parabolic problem. 

If we want to analyze now the case where the function is given by Ge{x) — a(x) + b{x)g{x/e) and the 
functions a, b are smooth but not necessarily constant, it is intuitively true that the limiting equation should 
behave like (1.3) with a diffusion coefficient depending on x somehow. Actually if we look at the thin 



domain in a small neighborhood of a point ^ G (0, 1), we will approximately see a domain with very high 
oscillations but locally the domain behaves like the thin domain with the function x — > a((^) + b{S,)g{x/e). 
Therefore, it is expected that if we freeze the coefficients of the limit equation in a fixed point ^ e (0, 1) we 
should recover equation ( 1.3 ). Although this intuitive argument will turn out to be true, it does not give us a 



complete information about the limit equation, specially, the way in which the dependence on x is explicitly 
stated in the limit equation. For instance, it is not clear at this stage whether the limit equation should be 
-{qix)wj:)x +w = / or -q{x)wxx + w = f with q{x) = \y*\x)\ ~ ^^^§^{yi,y2)}dyidy2, or maybe 

- \Y'\x)\ i'^i^)^^)^ + ^ f where r{x) = - ^^^^(yi, 2/2) 1^2/2, or maybe other. Observe that ah 

these equations coincide if wc consider the purely periodic case. 

In order to accomplish our goal and obtain the limit equation in the general case, we propose a technique 
that consists in solving first the piecewise periodic case, that is, the case where the functions a{x) and b{x) 
are piecewise constant and then do an approximation argument to obtain the limit equation in the general 
case. This is a subtle argument since we are approximating first the functions a and b by piecewise constant 



functions, say a (x) , b (x) so that 



a (a;) I + \ b{x) — b {x)\ < S and obtain the limit equation for (5 > 



fixed, passing to the limit as e — 0. Then, in this limit equation, which will depend on S, we pass to the limit 
as 5 — 0. But the limit we are interested in is taking first 5 — > for e > fixed, so we obtain the domain given 
by the function a{x) + b{x)g(x/e), and then we pass to the limit as e — ?■ 0. But, a priori there is no garantee 
that these two limits commute. We will actually show that these two limits commute by proving that the 
solutions of problem (1.1 ) in two domains fi'^ = {{x, y) : < x < 1,0 < y < Gi{x) = a{x) + b{x)g{x/e)} and 
51^ — {{x,y) : < X < 1,0 < y < Ge{x) = a{x) + b{x)g{x/e)}, are close in the norm uniformly in e when 
a, b and a, b are close. This result, which can be regarded as a domain perturbation result but uniformly 
in e, guarantee that the two limits commute and will show that the limit problem is given as above. We 
remark that this domain perturbation result is not deduced from standard and known results on domain 
perturbation techniques since we are able to compare the solutions of an elliptic problem in two families of 
domains which also depend on a parameter and the way this domains depend on e is not smooth at all. 

We strongly believe that this technique of solving first the piecewise periodic case and then use an 
approximation argument, uniform in the parameter e, can be used in other problems to obtain the appropriate 
homogenized limit for the non periodic case. 
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Following this agenda, we solve first the piecewise periodic case, that is, the case where the function 
Gi{x) = a^{x) + h^{x)g{x/e) and a'', are piecewise constant. We consider the points = < < • ■ • < 
£,N-i < = 1 E^nd assume that the functions and are constant in each of t he interval (Ci-iiCi)' ^^Y 
a'' (a;) = a^, b[x) = bi. We show that the limit equation is of the same form (1.3) in each of the intervals 
{£i~i,£i), that is, 

- QiWxx + w = f{x), a; e ,i^l,...,N (1.5) 

where 

It = I ^ (l - -g-^iyi,y2)^dyidy2, 



Y* 



and Xi is the unique solution of (1.4) in the cell Y* where 

= {(2/1,2/2) e : 0<7;i<L, 0<y2<a^ + b,g{yi)}. 



Moreover, equation ( 1.5 1 is suplemented with Neumann boundary conditions at a; = 0, a; = 1 and with some 
"matching" condition at the points ^i, i = 1, . . . , N — 1, which are continuity of the function and some kind 
of Kirchoff type conditions, see Section |3] for details. Actually, if we look at the variational formulation of 
the limit equation, we obtain 



' i^q^{x)p^{x) Wx{x)ipxix) +p^{x)w{x)ip{x)^dx = J p^ (x) f (x) ip dx 

where q^{x) = qi and p^{x) = in (6-1, Ci)- 

Now we will be able to pass to the limit in this equation as (5 — > and obtain the limit problem: 

1 



(1.6) 



p{x 



where 



■{r{x)wx)x + w = f{x), a;e(0, 1) 
u;'(0) = w'{l) = 
dX{x) 



(1.7) 



r{x) = p{x)q{x) = -r / |l ^-^(2/1, 2/2) 1^2/1^^2/2, 



\Y*{x)\ 
L 



p{x) = 



and X{x) is the unique solution of (1.4 1 in the basic cell Y*{x) which depends on the variable x and it is 
given by 

Y*{x) = {{yi,y2)eM? : < yi < L, < y2 < a{x) + b{x)g{yi)}. 



Equation (1.7 1 is the limit equation we were looking for. 



Finally, we would like to remark that although we will treat the Neumann boundary condition problem, 
we may also impose different conditions in the lateral boundaries of the thin domain i?"^, while preserving 



the Neumann type boundary condition in the upper and lower boundary. That is for problem ( 1.1 1 we may 
consider conditions of the Dirichlet type, w'^ = 0, or even Robin, + I3w'^ = in the lateral boundaries 
{(0,2/) : < y < eGe jO)} , {(1, y) : < y < eGe(l)}. The limit problem will preserve this boundary condition. 



That is, in problem (|1.7|) we will obtain the conditions w = or 4^ + f3w = at x = and a; = 1 



dN 



We describe the contents of the paper. In Section [2] we set up the notation, obtain some technical results 
that will be needed in the proofs and state the main result. In Section[3]we obtain the result for the piecewise 
periodic case. In Section [4] we show the continuous dependence result on the domain, uniform in e, that 
will be the key argument to obtain the limit problem. In Section [5] we provide a proof of the main result. 
Finally, we include an Appendix where we analyze the behavior of the basic function X solution of (1.4) as 
we perturb G. 
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2. Basic facts, notation and main result 

We consider the one parameter family of functions Ge : / — ?► (0, oo), where / = (0, 1), e £ (0, eo) for some 
eo > 0. We will assume the following hipothesis 

(H) There exist two positive constants Go, Gi such that 

0<Go<GAx) <Gi, Vxel, VeG(0,eo). (2.1) 
Moreover, the functions Ge(-) are of the type Ge(x) = G{x,x/e), where the function 

G:/xM ^(0,+oo) 

{x,y) ^G{x,y) ^''^^ 

is L-periodic in the second variable, that is, G(a;, y + L) = G(x, y) and piecewise with respect to 
the first variable, that is, there exists a finite number of points = < S.i < ■ ■ ■ < ^n-i < ^at = 1 
such that the function G : (^i, <^i+i) x M — (0, +cxi) is G^ and such that G, G^ and Gy are uniformly 
bounded in x M and have limits when we approach and 

One important example of a function satisfying the above conditions is 

N 

G{x, y) = a{x) + ^ br{x)gr{y) 

where the functions a, 6i,..,&7v are piecewise G^ in / = (0,1) and the functions 5i,..,(?Ar are all G^ and 
L-periodic. 

We define the thin domain as 

R' = {{xx,X2) G I e /, < 0:2 < eG,(a:i)}. 

In this work we study the asymptotic behavior of the solutions of the Neumann problem for the Laplace 
operator 

-Aw' + = in R' 
— — =0 on oR 

with /' e L'^{R'^) and where v'^ = {i^f, i^|) is the unit outward normal to and is the outside normal 
derivative. 

From Lax-Milgran Theorem, we have that problem (2.3) has a unique solution for each e > 0. We are 
interested here in analyzing the behavior of the solutions as e — >■ 0, that is, as the domain gets thinner and 
thinner although with a high oscillating boundary. 

To study the convergence of (2.3) in the thin oscillating domain R'^, we consider the equivalent linear 
elliptic problem 

d'^u' 1 d^u" „ . 

2 9 2+^ / mil 

dxi dx2 ^2 4) 

0x1 0x2 
where f € LP'iJY) satisfies 

WrWLHR^) < C (2.5) 
for some G > independent of e and now N'^ = (A^^, A^|) is the outward unit normal to dil"^ and il*^ C is 
given by 

= {(xi, xz) e I e /, 0<X2< G,{xi)}. (2.6) 
Observe that the equivalence between the problems (2.3) and ( |2.4[ ) is obtained by changing the scale of 
the domain R'^ through the transformation (x, y) — >■ (x, ey), (see JlO] for more details). Moreover, the domain 
is not "thin" anymore but presents very wild oscillations at the top boundary, although the presence of a 



6 



J. M. ARRIETA AND M.C.PEREIRA 



high diffusion coefhcient in front of the derivative with respect the second variable balance the effect of the 
wild oscillations. 

The domain il' is related to the ones analyzed in the papers [H [U |H] but the fact that in our case we have 
a very high diffusion in the y-direction makes our analysis and results different from these other papers. 

We write Hl{U) for the space H^{U) with the equivalent norm 



L\U) 



dw 
dxi 



+ 



1 



dw 

0X2 



LHU) 



given by the inner product 



A ^ f (a del) dip 1 d(j) dip I 



where U is an arbitrary open set of M'^, which may depend also on e. 
The variational formulation of (2.4) is find u'^ e H^{fl'^) such that 



ll^l^ + 4?^l^+^v}d^id^2= / r^dx,dx2,y^eH\n^). 

I oxi oxi 02:2 0x2 J Jsi« 



(2.7) 



Remark that the solutions u"^ satisfy an uniform a priori estimate on e. In fact, taking ip = u'^ in expression 

l^^lhin-) < WrhnnnhlLHan- (2-8) 
Wlhw) Ve>0. (2.9) 



(2.7), we obtain 



Consequently, 





1 


1 


du"^ 


dxi 




' e2 


dx2 



2 

L2(n=) 



du" 




du" 




and - 


< 


dxi 


L2(n«) e 


dx2 





Hence, it follows from (2.5) that there exists C > 0, independent of e > 0, such that 

II fii Ou'^ , 1 du"^ 

\ u \ L^iQ'), T, — and - - — 



< C. 



(2.10) 



Observe that the solutions tt *^ of (|2.4|) can also be characterized as the minimum of a functional. That is, 

(2.11) 



e H'^in") is the solution of if only if 

Vein') 



min VJip), 



where 



2 Jn' '-oxi 



dx2 



ip^^dxidx2 — J f'^ipdxidx2- 



It follows from (2.7 1 with ip ~ u'^ that 



Ve{u') = -\ I ru'dxidX2. 



Hence, due to (2.5) and (2.9) we obtain 



(2.12) 



Important tools for the analysis are appropriate extension operators for functions defined in the sets fi'^. 
We will obtain such operator and we will be able to construct it even for more general domains that the ones 
like ^\ 

Hence, let us consider the following open sets: 

O = {(xi,X2) e I a;i e / and < X2 < Gi} 
= {{xx,X2) e I e / and < X2 < G^ix^)] 
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where / C M is an open interval, : / M is a piecewise C^-function satisfying < Go < G^{xi) < Gi for 
aU a; e / and e > and there exists Q — < < . . . < ^n-i < (,n = ^ such that Gg is in the intervals 
(Ciifi+i)- Let us define 

77(e) = sup {\G:ix)\} (2.13) 

and assume r]{e) < +oo for fixed e, although in general 77(e) — > +00 as e — ^ 0. 
Also, we denote by 

d^O\ uiIi{(C„X2) : Go < X2 < Gi}. 

Notice that O" C O. 
Lemma 2.1. With the notation above, there exists an extension operator 

p, e c{LP{o'),LP{d)) n c{w^'PiO'),w^'P{d)) n ciWg'P{0'), Wg;P{d)) 

(where Wg'^^ is the set of functions in W^'^ which are zero in the lateral boundary di) and a constant K 
independent of e and p such that 



\\Peip\\Lp^O) ^ K\\(p\\lp^o') 



dPeV 



dxi 



LP{0) 



< K 



dip 



dxi 



vie) 



dip 



dPeif 



0X2 



LP{0) 



< K 



dip 



dxy 



dx2 

LP(0') 



LP{0')^ 



(2.14) 



(P^ip)(xi,X2) 



for all ip g W^^'^'(C^) where 1 < p < 00 and rj{e) is defined in (2.131. 

Proof Observe first that the set 0° = (0, 1) x (0,Go) C O" for all e. Hence, if we have that Gi < 2Go, 
which implies that Ge(xi) < 2Go, we can define the operator: 

p{xi,X2) {X1,X2) € O" 

pixi,2Ge{xi)-X2) {xi,X2) eO/0\ 

Observe that this operator is obtained through a "reflection" procedure in the X2 direction along the 
oscillating boundary. It is straigth forward to check that this operator satisfies (2.14 1. 

If we are in the case where Gi > 2Go, we will need to extend first the function p^ao in the direction of 
negative X2, with a finite number of successive reflections. That is, if ipo is defined in O"^ then we extend ipo 
to the set (0, 1) x (— Go,0) with the reflecting along the line X2 — 0. This produces the function 



ipi{xi,X2) 



(Pq{xi,X2) (a;i, 2:2) with < 2:2 < Gc(a;i) 
i^o(a;i, -a;2) (a;i, 2:2) with - Gq < 2:2 < 0. 



We can continue producing these reflections inductively as follows 

(Pn-iixi,X2) {xi,X2) with - {n-'l)Go < X2 <G^{xi) 

i^„_i(xi, -X2 - 2(n-l)Go) (a;i,2;2) with - nGp < 2:2 < -(n- l)Go 



ipn{Xi,X2)-' 



Choosing n large enough so that riGo > Gi, constructing (/?„ and applying to ipn the procedure by 
reflection in the 2:2 direction along the oscillating boundary, we obtain the extension operator which 
satisfies ( |2.14[ ). □ 

Remark 2.2. 1 ) This operator preserves periodicity in the Xi variable. That is, if the function p^ is periodic 
in xi, then the extended function P^Pe is also periodic in xi. 

2) This result can be applied to the case G^{x) — G{x) independent of e. In particular, we can apply the 
extension operator to the basic cell. 

Now we are in contidion to state our main result. We consider the general case, that is, the domain fi*^ is 
given as 



= {(2:1,2:2) e 



Xi 



el, <X2 <G,{xi)}. 
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where the function Ge(-) satisfies hypothesis (H). Recall that we denote hy fl — (0, 1) x (0, Gi) and 

n^n\ ufLim, X2) ■.Go<x2< Gi}, 

where the points Co: Cii • • ■ i are given by (H). 



Theorem 2.3. Let v!^ he the solution of (2.4-) with /"^ G L'^{Q,'^) satisfying ||/'^||l2(si=) < G for some positive 

constant C independent of e > 0. Assume that the function /'(a;) — J^'^^'^ f{x,y)dy satisfies that /'^ ^ /, 
w-L2(0,1). 

then 



2.1 



Then, there exists u G H^[^l), such that, if Pe is the extension operator constructed in Lemma 

P.u'^u w-H\n) 

where u{xi,X2) depends only on the first variable, that is, u{xi,X2) ~ u(xi), and u is the unique solution of 
the Neumann problem 



r{x) Ux{x) (fixix) +pix) u{x) (p{x)'jdx = J f{x) (p{x) 



dx 



for all (fi e H^(I), where 



r(x) 



dX{x) 



idy2 



(2.15) 



(2.16) 



and X(x) is the unique solution of problem 



L 

-AX(x) = m Y*{x) 
dX{x) 



dN 
dXjx) 
dN 



on B2{x) 



Ni on Bi(x) 



(2.17) 



X{x){0,y2) = X{x){L,y2) on Bo{x) 
X{x) dyidy2 = 

Y'(x) 

in the representative cell Y*{x) given by 

Y*{x) = {{yi,y2) e : < 2/1 < L, < < G{x,yi)}. 



(2.18) 



Bo{x) is the lateral boundary, Bi{x) is the upper boundary and B2{x) is the lower boundary of dY*{x) for 
all X e /. 



Remark 2.4. i) If the funct ion r (x) is a continuous function, then, the integral formulation {2.15) is the 
weak formulation of problem (1.1) with f{x) = f{x)/p{x) 

ii) If initially the function f'^[x,y) = fa{x), then it is not difficult to see that f ix) = Ge{x)fo{x) and 



G,{x) 



\Y'{x)\ 



p{x) as e — !■ 0, and therefore, f{x) = p{x)fo{x). 



3. The piecewise periodic case 



In this section we find the limit of the sequence {u'^}e>o given by the Neumann problem (2.4) as e goes 
to zero for the case where the oscillating boundary is piecewise periodic. 

So let us consider that the family of domains Jl^ satisfies (H) and morever the function G is independent of 
the first variable in each of the domains (Ci-i, Ci) 'x^- That is, there exist = < < • • ■ < ^n-i < £,n = ^ 
so that the function G from (2.2| satisfies that G{x,y) = Gi{y) for x E Ii = and Gi{y + L) = Gi{y) 
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for all y E M.. Moreover, the function Gi{-) is for alH = 1, 2, . . . N, there exists < Gq < Gi such that 
< Go < Gi(-) < Gi for all i = 1, . . . , A^; and the domain is given by 



|(a;,y) : < x < ^,,0 < y < G.,{x/e),i = 1,.. .,iv| 
U^7'{(e.,y),0 < y < min{G,_i(6A),G,(6A)}} 



U 



(3.1) 



(see Figure [2]). Denote also by fl the open rectangle = / x (0, Gi) and by 

n = n\ uili{(e., X2) ■.Go<x2< Gi}. 

Observe that for Q'^ and f2 we have the extension operator constructed in Lemma 



2.1 




Figure 2. Piecewise periodic domain ff^ 



We can show, 



Theorem 3.1. Assume that e L'^{Q'^) satisfies (2.5) and the function f'^{x) = /q'^'^^' f{x,y)dy satisfies 
that /"^ ^ /, w-L^(0, 1). Let be the unique solution of (2.4-)- Then, there exists u G H^(il) such that if 
is the extension operator constructed in Lemma \2.1\ we have 

P^u'^ii w- H^{n), s- L'^iCl) 

where u{xi,X2) = u(xi) in Q and u{-) is the unique weak solution of the Neumann problem 

for all ip e H^{I), where p{x) and r{x) are piecewise constant functions defined as follows: p{x) — pi for all 
X e (6-1, Ci) where 

\y; 



L ' 



i = 1. 



.N 



where Y* is the basic cell for x G (Ci-iiCi); that is 

Y* = {(2/1,2/2) e M2 : < yi < L, < j/2 < G,(yi)} 
and r{x) = ri for all x £ where 

dX, 



(3.3) 
(3.4) 
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and the function Xi is the unique solution of 



-AX, = m Y* 



dX, 
dN 



= A^i on B\ 



(3.5) 



X,{Q,y2) = X,{L,y2) on B}, 
Xi dyidy2 = 



(3.6) 
(3.7) 



Remark 3.2. If we define /o(a;) = f{x)/p{x), then problem (3.2) is equivalent to the following: 

~ qiUxx{x) +u{x) = /o(a;) x G (Cj-i,Ci) 
for i — 1, ---tN, where qi — ri/pi, satisfying the following boundary conditions 

[ QiUxiCi-) ~ qi+i Ux{ii+) = i = 1, ....,N - 1. 
Here, Ux{£,i'^) denote the right (left) -hand side limits of at ^i. 

Proof. Let us consider the family of representative cell F/, i = 1, 2 . . . , TV, defined by 

Y* = {(yi,y2) e : < yi < L and < < G,{yi)} 

and let Xi be their characteristic function. We extend each Xi periodically on the variable yi € M and denote 
this extension again by Xi, for « = 1, . ■ • , -/V. 

If we denote by xl tbe characteristic function of the set f)- = {{x,y) : < a; < ^i, < y < Gi{x/e) — 
G(a;,x/e)}, we easily see that 



X\ 

x'i{xi,X2) ^Xi{ — ,X2), for (xi,X2) e ^l- 

Let us also denote by fi* the rectangle 17* — {(x, y) : < a; < ^i, < j/ < Gi}, for i = 1, 
Let us also consider the following families of isomorphisms : A], i— >■ Y given by 

Tk(Xl,X2) = ( ,X2j 

where 

Al = {ixi,X2) e I ekL < xi < eL{k + 1) and < a;2 < Gi} 
r = (0,L) X (0,Gi) 

with k eN. 

Let us consider the following auxiliary problem given by 

-AX, = in Y,* 
dX, 



(3.8) 



(3.9) 



dX, 
dN 



= on Bh 
ON ^ 

G',{yi) 



on B\ 



X,{Q,y2)=X,{L,y2) on 



(3.10) 



Xi dyidy2 = 



Y* 



where Bq is the lateral boundary, Bl is the upper boundary and B2 is the lower boundary of dY* . 
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Taking the isomorphism (3.9 1 and the family of extension operators 

eC{H\Yn,H\Y))n£ [L^ {Y* ) , (Y) ) 



defined by Lemma 2.1 with Gc{xi) — Gi{xi) independent of e and Y = (0, L) x (0, Gi), see Remark 2.2 
define the function 



we 



Xi-e{P'X,i 



,xi — eLk 



,X2)), for {xi,x2) e n'nAl, 



I = 1, 



,7V. 



Notice that this function is defined in uflj^Jl' and it is well defined. For e > fixed and for (xi,X2) € fi* 
for some i = 1,. .. ,N then there exists a unique k G N such that {xi,X2) £ A'j^. Observe that uflj^O' = 
^l\U^=i^{{^i,y) : < y < Gi} and that e H'^{ufLi^")- 

We introduce now the vector rj'^ — [rjl, 77I) defined by 



Since gf^ 



e c3i;i c3a;2 



?7^(a;i,a;2) = —{xi,X2), {xi,X2) e U.^^rj*, r 

= -rf- we have that 
oy2 



1,2 



dXi xi~ ekL dX Xi 

Vi[xi,x2) = 1 - ^— ( ,2:2) = 1 - ^(—,2:2) := m{yi,y2), 



, , dXi xi — ekL dX xi 

'q^(xi,X2) = -£75—1 ,X2) = -e7^( — ,a;2) 

01/2 e o?;2 e 

for (yi,y2) = (^^^,2:2) e K,* and ixi,X2) € for i = 1, . . . , TV. 
It follows from definition of X that the functions rjl and 77^ satisfy 

|^ + 4|^-0inl^^z = l,...,iV 
axi e-^ 0x2 



r,lN', + ^v'^NI = on (xi,G,(^)), * - 1, . . . ,iV 



'yj^i' + ^^2^^2=0 on (a;i,0) 



where 



N' = iNl,N^) 



TV" = (0,-1) on (a;i,0). 



In fact, by (3.10) we have 
dxi 3x2 



l(d^X, xi d^X, xi \ „ . . . 



(3.11) 



(3.12) 



(xi,G,(^)),* = l,...,iV 



GK^)(l-||(^.^2))+||(^,X2) 

ye^+y(W 

G^(^) + (i7('r-^2), li(^,2;2)) • (-GU^^), 1) 



G' 
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on (a;i,G,(3^); 



for i = 1, 



, TV. Moreover, 



r^lNl + \rjlN^, = (1 - • + -^(^,^2) = on (:ri,0). 



Therefore, multiplying (3.12 1 by a test function 1^ G H^{VL) with (/? = in neighborhood of the set 
UiLo{(Cii ^2) : < X2 < Gi} and integrating by parts, we obtain 

(dill 1 ar^n 

((5 7^ h — dxidx2 

\dxi e^dx2j 



^ fyyJTV^ + ^il'^Ni] dS ^ f ( + 



dxidx2 



That is 



, 9V9 1 ^(^ 



(3.13) 



Let ( 
obtain 



.Ti) S C^(U^Q^(^i, ^i+i)) and considering the test function ip = ip lo'^ m. (2.7) and in (3.13), we 



f'^{4> uj'^)dxidx2 



{ dxi dxi ^ 
r du'^ d 
I dxi dxi 
d 



1 du' d 
8x2 dx2 
1 du" d ^ 
8x2 8x2 
8 



-Vi 



uj'^) + u'^{4> w*^) I da:; 1 da; 2 

w^) + ^^^^^('^ + ^"^^ w'')|cia;ida;2 

(8 18 ~i 

V^dxi^^ + ^V2-Q^{(l^^ljdxidx2 

( 8u'^ 8(j) , 8u^ du'^ , 1 8u'^ 8uj^ , , , , 

!> 1 1 

- ^it; — ^^2^ — ^da::ida;2. 

1 8x1 0x2 J 



8x 



Using that rjl — we cancel the appropriate terms and obtain 



I 8x1 8x1 



Vi-K—u 
0x1 



(f) u:'^^dxidx2 

r<i>oj'dxidx2, V0 e Co"(u^oite,e.+i)). 



(3.14) 



(3.15) 



On the other hand, we have obtained before the weak formulation of problem ( |2.4[ ), that is 

/ ll^l^ + 4l^l^ + "''^V^i^^2= / rcpdx,dx2, y^p e H\n^). (3.16) 

Jo= 0x1 8x1 da;2 da;2 J Jo-^ 

Now we need to pass to the limit in (3.15) and ( |3.16 ). In order to accomplish this we need to write both 
expressions as integrals in the same domain. For this, we will use the extension constructed in Lemma 
2.1 the standard extension by zero, that we denote by~ and the characteristic function of il*^ as follows: 
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Jn 

I \^^+x'P.u'^]dx^dx2^ [ xT'pdxidx^, y^eH^O,!). 
Jq I dxi dxi J 



(3.17) 



(3.18) 



Observe that in this last equahty we have taken (p € H^{0, 1) and the term inchiding partial derivatives with 
respect to X2 do not appear. 

We want to pass to the limit in the expressions above, (3.17) and (3.18) . In order to accomplish this, we 
pass to the limit in the different functions that form the integrands. 

(a). Limit of Xe- 



From (3.8), we have for i = I, . . . , N , 

Xt{-,X2) e,{x2) y 



X^is,X2)ds u;*-L°°((e._i,6), V.T2e(0,Gi). 



(3.19) 



Observe that the limit 9i does not dependent on the variable xi € (^i_i,^i), although it depends on 
i = 1, . . . ,N. Moreover, we can get the area of the open set Y* with the formula 

r-Gi 

e,{x2)dx2 = \Y*\. (3.20) 



L 



Also, from (3.19) we have that 



(p{xi,X2) |xi(a;i,a;2) - 9iix2)^ dxi ^ as e ^ 

1 

a.e. X2 e (0, Gi) and for all ip £ L^{n). So, due to 

ip{xi,X2) ^Xi{xi,X2) ~ 0i{x2)^ dxidx2 = / H-{x2)dX2 



and \H-{x2)\ < / \(p{xi, X2)\dxi, 
we can get by Lebesgue's Dominated Convergence Theorem that 

where 9{xi,X2) = Oi{x2) for xi € (Ci-i,Ci), i = 1, 2, . . . , TV. 
(b). Limit in the tilde functions 



(3.21) 



Since is uniformly bounded, we get from (2.8 1 that there exists M independent of e such that 



1'"11l2(o), 









and - 


dxi 


L2(o) e 


0X2 



L2(a) 



< M for all e > 0. 



(3.22) 



Then, we can extract a subsequence, still denoted by m^, and such that 

2/- 



dxi 



C w-L^{n) and 



(3.23) 



dxo 



s-L^{n) 
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as e — ^ for some u* and ^* e 

Moreover, since ||/'||l2(q,-) < C independent of e, we have ||/^||L2(f2) < C and therefore, the function 
defined by 

/•Gi 

(3.24) 



Jo 

satisfies that /'^ G L^{0, !)• Hence, via subsequences, we have the existence of a function / = f{xi) E L^{Q, 1) 
such that 



-i'(o,i). 



(3.25) 



(c). Limit in the extended fun ctions 

Using the a priori estimate (2.81, the fact that G H^(il'^) and using the results from Lemma 
extension operator we get that 



2.1 



on the 



\P.u' 



dP 



dxi 



and - 

L2(o) e 



dP,; 



dX2 



< M for all e > 



where M is a positive constant independent of e given by estimate (3.22) and Lemma 
extract a subsequence, still denoted by P^u'^ and a function uq S H^{n), such that 

Peu'^uo s-L'^{n) 
and^^^O s-L2(m. 

OX2 



2.1 



(3.26) 

Then, we can 

(3.27) 



A consequence of the limits (3.27) is that uo{xi,X2) does not depend on the variable X2- More precisely, 

duo 



dx- 



■{xi, X2) — a.e. f2. 



(3.28) 



In fact, for i = \, . . . ,N and for all ip e C^i^i), we have by (3.271 that 

dip 
dx2 
dPM" 



uo dxidx2 



lim / P.u" 

e-i-O 



lim 

€-»-0 



8X2 



dxidx2 
(p dxidx2 = 



which implies that uq{xi,X2) does not depend on X2- Morever, since the rectangle / x (0,6*0) C il*^ for all e 
and e H^{Ix (0, Go)) we have from ^3.27^ that uq w - HHIx ( Gp) ) and therefore uq G iJ^O, 1). 

Also, we note that = x'^PeU'^ a.e. fl. Thus, it follows from (3.21), (3.23) and (3.27) that we have the 
following relationship between u* and Uq 



u*{xi,X2) — 9i{x2)uo{xi) a.e. {xi,X2)&^i, i — l,...,N. 



(3.29) 



(d). Limit in cj^. 

With the definition of cj^, we have for alH = 1, 



jw*^ — Xi\^dxidx2 — 



Ainfi' 



Yi 



'\{PX,){y,,y2)\^dy,dy2< / Ge3|X,(yi, 2/2)|'d2/idy2 



and so, 



/ Iw' -a;ipdxidx2 « V / Ce^\X,{yi,y2)\^dyidy2 
w J C\X^{yi,y2)\^dyidy2 ^ as e ^ 0. 
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Similarly, 



dxi 



(w" - Xi) 



dxidx2 



d{PXi) 



< e C 



dyi 



and 



(w" - Xi) 



dxidx2 



Yi 



dyi 

dy2 
dX, 



{yi.yi) 



c 



dy2 



iyi,y2) 



{yi,y2, 



edyidy2 



dyidy2 



edyidy2 



dyidy2 



Also, we have 



f r f OX' 

J dx^^^"'^^^ dxidx2f=i^^j Ce\-^{yi,y2)\^dyidy2 



C 



dX, 
dyi 



(yi,2/2) dyidy2 



for all e > and 



Hi 



dx' 



(w' - Xi) 



dxidx2 < e / C 



dy2 



(yi,2/2) 



dyidy2 — J' as e — >■ 0. 



Then, we can conclude 



and 



uj'^xi s - L'^in) and w - H^{n,), i = l,...,N, 



dx-. 



->0 s~L\n). 



(3.30) 



(e). Limit of rj\ 

Let rj'^ = ry'x'^ be the extension by zero of the vector -q'^ to the whole Vl. We can obtain by the Average 
Theorem that 

'-{s,X2)]xi{s,X2)ds ri{x2) - Ci) (3.31) 



Vl{xi,X2) 



1 



L 



1 - 



dyi 



where Xi is the characteristic function of Y* . 
Hence, arguing as (3.211 we can prove 

Vl-^f w*-L°°(n). 
where ■r{xi,X2) = fi{x2), {xi,X2) £ fli, i — 1, . . . , N. 



(3.32) 



Now, by the convergences shown in (a)-(e) above, we can pass to the limit in (3.17) and in (3.18). We 
obtain, for all G Co°°(U^-\e», 6+i)), 



j {^*^^ -^1 ~ "0 + ^"o0a:i|(ia;i(ia;2 = y f{xi)(l){xi) xidxi. 
Observe that gf-l'/'^^i) = ^*^^'§T' + ^*'^- Consequently, we have 



|r^^('/'2;i) - 0r - f-Q^ uq + 9uo(f)Xi^dxidx2 



f{xi)(j){xi) xidxi 



(3.33) 
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for all cf> e C5-(uto^(e.,6+i)). From (|3J8f, we get 

dip 



(3.34) 



In particular, via iterated integration and (3.34 1, we get 



d^ixi) , \Y*\ 



dxi L 

fixiMxi)dxi, y^eH\o,i). 



uo{xi)ip{xi)^dxi 



(3.35) 



Taking (p — cjjxi in (3.34), we get 



|C*^-(0a;i) + 6'wo0xi|da;ida;2 = J f4>xidxi. 
Hence, it follows from ( |3.33[ ) and ( |3.36[ ) that, for all (f> e C^{^^=o^{S.^,^^+l)), 

= j ^4>C +f^^UQ^dxidx2^ J ^(j)^ - r(f)^^^dxidx2 

where we have performed an integration by parts to obtain the last integral. Observe that this integration 
by parts can be performed since e C^(U^Q^(^i, ^i+i)) and f does not depend on xi in each of the domains 
rii. Hence, if we define 



(3.36) 



(3.37) 



f Iff dX- "I 

j fi{s)ds=-J i^l--^{yi,y2)jdyidy2, i = l,...,- 

r{xi) = ri, for a;i e (^i_i,^i), i = l,...,N 
re get 

)( y^"^f (xi,X2)dx2-r(xi)^^^)da;i=0, V0 G Co-(u£oi(e„ Cm)) 



and we denote by 

and performing an iterated integration in (3.37) we get 

r-Gi 

(Xl 



which implies that 



C{xi,X2)dx2 = r(xi) ^"°^^^\ 
ox I 



a.e. Xl e (0, 1). 



Plugging this last equality in (3.35) we get 

\Y* 



N 



' duQ dip 
^ ' dxi dxi 



uq ipdxi 



f^dxi, V^€i7i(0,l). 



(3.38) 

(3.39) 

□ 



4. A DOMAIN DEPENDENCE RESULT 



In this section we are going to analyze how the solutions of (2.4) depend on the domain i}'^ and more 
exactly on the function G^. As a matter of fact we will show a continuous dependence result with respect 
to the functions G^. 

More precisely, assume and Ge are piecewise continuous functions satisfying (2.1) and consider the 
associated oscillating domains 17*^ and f^*^ given by 

= {{xuX2) e I Xl e /, 0<X2< G,{xi)} 

Cl" = {ixi,X2) e I Xl e /, < X2 < G'e(xi)}. 



Let u*^ and u'^ be the solutions of the problem (2.4) in the oscillating domains fl"^ and f2'^ respectively with 
e L^(M^). Then we have the following result: 
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Theorem 4.1. There exists a positive real Junction p : [0, oo) t-^ [0, oo) such that 

with p{5) — > as (5 — > uniformly for all 

• e > 0; 

• piecewise functions and with < Go < G^{-),G^{-) < Gi and 

llG.-G,|Uoo(o,i) <<5; 

Remark 4.2. The important part of this result is that the function p{5) does not depend on e. Only depends 
on Go and Gi . 

To prove this theorem, we use the fact that u'^ and u'^ are minimizers of the quadratic forms 

1 /" d ^ 1 d ^ P 



^'^^^^2 LAd^ +^'}d^^dx,- f^^dx.dx,. 



(4.2) 



That is, we have 



VJu') = min VJifi) 
= min V^{ip). 



In order to prove Theorem 4.1 we wiU need to consider the minimizer of the functionals Vg, Vg, and plug 
them in the other functional. For this, we need to transform the function u'^ into a function defined in fi*^ 
and the function m*^ into a function defined in Sl*^. One possibility is to use some kind of extension operator 
as the one we have constructed in Section [2j But the problem with this approach is that the norm of the 
extension operators will depend on the derivatives of the function G^ and Ge and therefore it will be very 
unlikely to prove a results that will depend only on the L°° norm of G^ — G^. 

In order to transform the function u'^ (resp. u'^) into a function defined in fi*^ (resp. fi*^), we construct the 
following operators: 

Pi+r,:H\U)^ H\U{l + ri)) 

( X2 \ (4.3) 

{Pl+rj<P){xi,X2)='f\Xi,^-^j {xi,X2)eU 

where 

U{1 + 77) = {{xi, (1 + rj)x2) e I (a;^, a;^) e U} 

and ?7 C is an arbitrary open set. 

We also consider the following norm in H^{U) 

and we can easily see that 

MlHliU) = \\Pi+vf\\Hl,^^{U{l+r^)) (4-5) 

and 

Y^Mhiu) ^ Mhi,^,{u) < (1 + v)\Mh(uy (4-6) 

We have the following preliminary result about the behavior of the solutions near of the oscillating 
boundary. 
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Lemma 4.3. Let be the solution of the problem (2.4-) and let Pi+,, be the operator given by (4.3). 
Then exists a positive constant C = C{Gi, ||/'^||l2) independent of e £ (0, 1) such that 

ll^^llffiCO^Vn'CjJ^)) + ll-f'l+'j"1lHi(0'(l+r,)\f2») + II"' ^ -fi+fj"1li/i(o-) - ^Vv 

for all rj > 0. 

Proof. Since ry > 0, we have 17^ (y^) C ^l^ So, we obtain, 

1„ 



Vein') 



1, 



ll"'ll//i(n=\n^(^)) + 2II"'IIhH^2=(t^)) " L ^^i'^^' 



2 1 

lfl-i(f2'\f2^(TT^)) ^ 2 



" /^^ f''u'dxidx2 



> -\\u 



ll"'lll/i(n=\n^(^)) + ^^^^^WPi+v^'Whi^') 



f'^u'^ dxidx2 



where we have used (4.5), (4.6 1. Moreover, 

\\Pi+v^'\\Hiin^) = ll-Pi+^"' - + "1lffi(n« 



-Pi+>7"1Ih1(0= 

= ll-Pi+j)"' -"1lHi(n«) + ll"1lHi(n=) +2(^1+7,"' 

||Pi+„u^ - u^ll^wn.) + Il"'llffi(f2.) + 2 / (Pi+^u' - dxidxz 

satisfies the variational formulation (2.7) with ip — Pi^^u'^ — u'^ e H^{^1'^). 
rom (1471) that 



where we have used that u'^ satisfies the variational formulation (2.7) with ip — Pi+n'* 
Consequently, it follows from (4.7) that 

> i|l7/,^l|2 + IIP, .,,,1^ 



2ll"'llHi(f2^\o«(^)) + 2(1 + 77)"^'+""' " "'ll?^cMn=) 
1 /■ / 1 



/"^Ue dXidX2 



Hence, we obtain 



^ 2' 

+ -^V,{u')+ f (-^Pi+^u'-uA rdxidX2. 



2ll"'llHi(f2=\n = (T^)) + 2(l + 77)"'^'+''"' " "iHi(f2') 



< — — 

~ l + r] 

analyze the integral 

" ' 1 



V 



Now, we 



■^("') + / (t^Pi+vu' ~ u') r dxidX2. 



[ [ -^—Pi+T -uA f''dxidx2. 



To this, observe that 

M^(a;i,X2) - {Pl+rjU''){xi,X2) = u''{xi,X2) - u" (xi, = j ^ 
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which imphes 



u'^{xi,X2) — Xi, 



X2 
1 +?7 



< 



< 



ds 
ds 



{xi,s) 



ds 



1/2 



(xi,s) 



ds 



VX2 



1/2 



Thus 



u''{xi,X2) - u" Xi, 



_X2_ 
1+7? 



dX2 



< 



(xi,s) 



and we have 



Consequently 



llu' - Pl+,,w'^||L2(f2.) < 



1 + ?; 



dxo 



Pl+r^u"^ — I /"^ dxidx2 



1/2 



Gi. 



< 



< 



/ l/^M^I dxidx2 + f \{Pi+^u' - u') r\ dxidX2 



1/2 



dx2 



Gi. 



Then, it follows from (4.8) that 

1. 



£||2 



< 



1 /2 



Gi. 



Hence, due to (2.9), (2.12| and (4.10), we obtain 



\Hi(n'\n'(^)) 



|lPi+,X-"'lllfi(n. 



On the other hand, we have 
1, 



/'u' dxidx2 



2 ir'^^llffKo-i) ' 
2ll-Pi+';^1lHi_i+^(n=(i+^)) 



/"^u*^ dxidx2 



1, 



^'i+»7^1Ihi^^^(o^(i+„)\o=) + 2ll^i+'?"'llff.\i+,(o-) 



/'^m'^ dxidx2 



> 



(||i'i+,,u'liffi(n.(i+^)\n.) + ll-Pi+,,ulffi(o.)) - / fu'dxidx2 



2(1 + 7?) 



(4.9) 



(4.10) 



(4.11) 
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So, we can proceed as in (4.11) to get 



Putting together (4.11 1 and (4.12), we complete the proof of the lemma. 



We are in conditions now to proof the main result of this section. 



Proof of Theorem ^.1. If we define 77 = 5 /Go, then under the hypotheses on and 
G(\\l°°(q.i) < S we have that 



Applying (4.13) and Lemma 4.3 we easily get that 



e||2 



So we will concentrate in the first term of (4.1). Therefore, 



= \\\Pi+vu'\\Hi{n^)- f'{Pi+^u'')dxidx2 

f'^Pl+rjU'^ dxidx2 
/ f'^Pl+rjU'^ dxidX2- 



But from (4.51 and (4.6 1 we get 



Moreover 



/ Pij^y^ii'^ dxidx2 — I f'^{Pij^y^u'^ — ii^) dxidx2 + I f^v!^ dxidx2 
JCi' JCi' Jh' 

but from Lemma |4.3t 

I / f\P^+.,u^-u,)dx^dx2\ < ||riL.(fV)ll^i+.'i' < Cri'"- 

Also, 

I / f^P,+r,u'dx,dx2\ < \\flmfiJPi+v^'\\mn,\n^) < ^V^'^ 



where we have used Lemma 4.3 and (4.13) 



Hence, putting all this information together, we get 

KK) < (1 + r,)V,{u') + iim\L^n^)\mmh') + Cr^'" < t{uc) + Cr^''\ 

Therefore 



On the other hand, we obtain by (4.4), (4.5) and (4.13) 

1„ 



/^u'^ dxidx2 
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But, due to (4.9) and Lemma 

1 



- 2(1 + 7,)' 

+2(Pi+,,u" - '"',"')iji(o«)) ~ y fu^dxidx^ 
+2 / {Pij^y^u'^ ~ u'^)f'^ dxidx2\ — I f^u"^ dxidx2 

- 2(rT^IIW-^i?.,(,, + ^v;(.^) 

Pi+nu" f dxidx2 - / fu" 



1 



1 + V Av 
4.3l we have 



^■^-"^ dxidx2- 



< 



< 
< 
< 
< 



1 



dxidx2 — / f'^u'^ dxidx2 



^ dxidx2 



1 


+ 


V 




1 




1 


+ 


V 




1 




1 


+ 


V 




1 




1 


+ 


V 




1 





■ / Pi+r^U^ dxidx2 — i f^u^ dxidx2 



i')\\r\\L2(n'{i+'n)\n') + Cqt]^^'^^ 

lendent of f a,nd n. 



that 



with Co, Ci and C2 independent of e and 77. 
Consequently 

Thus, it foUows from ( |4T6| and ( [ils] ) 
which imphes 

\\P,+,u^-u^\\l,^^^^<C^5''' 

J mdenenrlent ot f and 77. 

But, 



with C3 independent of e and 7;. 

But, from Lemma 4.3 we have that \\u'^ — Pi+^Ti^H^ij-jj,^ < Cy/fj. Hence, from (4.19) we get 



(4.17) 



(4.18) 



(4.19) 



(4.20) 
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for C independent of e. This provers the result. q 



5. The General Case 
We provide now a proof of the main result, Theorem |2.3| 



Proof. From estimate (|2.10 ) and (|2.14|) we derive 

L2 (o) 



\P.u' 



dP.u" 




dP^u^ 


and - 


dxi 


dx2 


L2(o) e 



L2(n) 



< M for all e > 



where M is a positive constant independent of e. Hence, there exists uq G H^{(1) and a subsequence, still 
denoted by PeU'^, such that 

Peu'^uo s-L'^{n) 
dP.u" 



(5.1) 



and 



dxo 







It follows from (5.1) that uq{xi,X2) = Uo{xi) on Cl. 



We will show that uo satisfies the Neumann problem (2.15). To do this, we use a discretization argument. 

Let us fix a parameter S > small enough and consider a function G^{x,y) with the property that 
< G^{x, y) — G{x, y) < S for all (x, y) € / x E and such that the function satisfies hypothesis (H) and 
is piecewise periodic as described in Section [3j To construct this function, we may proceed as follows. The 
function G is uniformly in each of the domains x ^ and it is also periodic in the second variable. 

In particular, for S > small enough and for a fixed z e (Ci-i-^i) we have that there exists a small interval 
{z — ri,z + 1]) with 7] depending only on S such that \G{x,y) — G(z,y)\ + \dyG{x,y) — dyG{z,y)\ < S/2 for 
all a; e {z — rj, z + ?/) n and for all y e M. This allows us to select a finite number of points: = 

< < . . . < <e[_i = 6 such that Cr_i - < V and therefore, defining G^{x, y) = G(Cr_i, y) + 6/2 
for ah X e (Cr-i,^!^!) we have that < G\x,y) - G{x,y) < 6 and \dyG^x,y) - dyG{x,y)\ < d for all 
(a;, y) £ ^i) x M. This construction can be done for alH = 1, . . . , iV. 

In particular, if we rename all the points S,^ constructed above by = zq < zi < . . . < = 1 (and 
observe that m — m{5)), the function G^{x,y) — Gf{y) with (x, y) £ {zi-i,Zi) x 
and L-periodic. 

We denote by Gf(a;) = G^[x,x/e) and consider the domains 

VL'^-^ ^{[x,y)£M?\x£l, 0<y<Gf(x)} 



i = 1, 



and G? is 



/X (0,Gi) 



Since Gf satisfies the hyphoteses of the Lemma 



= n'\ u™TH(^«,y) e | Go < y < Gi}. 

there exists an extension operator 



2.1 



1/6. 



satisfying the uniform estimate (2.14) with ry(e) 

Since £ ^^(rJe) with ||/'^||L2(f2«) < G, then if we extend this function by outside il'^ and denoting 

the extended function again by and using that Gs > G, we have that fg{x) = j^"^^^ f'^(x,y)dy = 
jGt(a;) jt^^^^y^^y _ /'(x) and by hypothesis, we have that fl = f^-^f w-L^(0, 1) 



solutions u"'^ of (2.4) in Jl''*^ satisfy 



Therefore, it f ollow s from Theorem 3.1 that for each 5 > Q fixed, there exist £ H^{0, 1) such that the 



(5.2) 



where u^{xi,X2) = u^{xi) on is the unique solution of the Neumann problem 
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{/(x) It* (x) ip^ix)+p^{x) u^ix) ip{x)^dx fix) tp{x) dx, V<y9 e 

where : 1 1-^ R and : / M are strictly positive functions, locally constant, given by 

{r^{x) = ri^s{x) 

where 



(5.3) 



(5.4) 



(5.5) 



The function Xi is the unique solution of (3.5 ) in the representative cells Y* defined in (3.4) for alH = 1, to. 
Now, we pass to the limit in (5.3) as 5 — » 0. To do this, we consider the functions and defined in 



X e / and the functions r and p defined in (2.16). We have that and p^ c onverge to r and p uniformly 
in /. The uniform convergence of to r in / it follows from Proposition |A.l| proved in the Appendix. The 
uniform convergence of p^ to p follows from the uniform convergence of G*to G as 5 ^ 0. 

Since we have the uniform convergence of and p^ to r and p respectively, we have by [H p. 8] or fS', p. 
1] the following limit variational formulation: to find u E H^{I) such that 



(5.6) 



' ^r(x) Ux{x) ipx{x) + p{x) u{x) ip{x)'^dx = J f{x)ipdx 



for all ip e H^{I). 

Hence, there exists u* £ H^{I) such that 



u* in H\l) 



(5.7) 



where u* is the unique solution of the Neumann problem (5.6) 
To finish the proof, we need to show that u 
Let us consider the open square flo 



uq in /, where uq is the function obtained in (5.1 ). 
/ X (0, Go) which satisfies Uq C fig for all i5 and e small enough. 



Observe that \\u* — ito||^2(/) = Gq ^\\u* — ""01112(5^^^) and therefore, to show that u* 



Uq it is enough to show 

that \\u* — '«o||'^2(Qp) = 0. Hence, adding and substracting the appropriate functions and with the triangular 
inequality, 



\\u* 



\\u* -u'\\L2^n„) + \\u' ~u^''\\mno) 



(5.8) 



for all e and S > 0. 

Now, let 7y be a positive small number. From (5.7) and from Theorem 4.1 we can choose a 6 > fixed and 
small such that — u^\\L^(no) ^ V a^id — u'^\\L^{Qf,) < rj uniformly for all e > 0. Now, from (5.2 1 for 
this particular value of S, we can choose ei > small enough such that — u'^'^\\L^{no) < for < e < ei. 
Moreover, from (5.11 we have that there exists £2 > such that — uo||l2(j7jj) < for all < e < £2- 
Hence with £ — min{£i,£2} applied to (5.8), we get — 1*011^2(^2^) < 4r]. Since rj is arbitrarilly small, then 
u* = Uq. □ 



Appendix A. A perturbation result in the basic cell. 
In the proof of the main result in Section [s] we have used the convergence of 



— > r, where r — 



n,5 



in the interval (zi_i,Zi), and ,5, r are defined by (5.5) and (2.16), respectively. In order to proof such a 
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convergence we need to analyze how the function X , solution of 

-AX in Y* 



ax 



5X 



on B2 
d2G 



on Bi 



on the representative cell 



1 + id2Gy 
= X{L,y2) on Bq 

X dyidy2 ^ 



Y* - {{yi,y2) eR^\0<yi<L, < y2 < G(yi)}, 



(A.l) 



depends on the function G. 

We will consider the following class of admissible functions 

A{M) = {Ge Ci(E),L- periodic, < Gq < G(-) < Gi,\G'{s)\ < M} 



(A.2) 



(A.3) 



and we will denote by Y*{G) and X{G) the basic cell (A.2) and the solution of (A.l) for a particular 
G e A{M). Observe that for each G € A{M), we have the extension operator Eq : H^{Y*{G)) -J> H^{Y) 
as constructed in Lemma 



2.1 



which will satisfy WEguWh^iy) < G\\u\\hi{y-'{g)) with G = C{M), but G 
independent of G, and therefore this constant can be chosen the same for all G G A(M). 

For each G £ A(M), we can consider the basic cell Y*{G) defined by (A.2). Here, we proceed as [21 p. 
84] and [10] . we begin by making the following transformation on the domain Y*{G) 



where F(z) 



L(5 : Y*iG)^Y*{G) 

(21,22) {zi,F{zi) Z2) = (yi,?/2) 
^ft^. The Jacobian matrix for this transformation is 

G(z) 



JLa(zi,Z2) 



1 

F'(zi)z2 F(zi) 



Gjz) 
G(z)- 



and observe that its determinant is given by \ JLq{zi, Z2)\ — F{zi) 
Using Lg, we can show that problem ( A.l ) in Y*{G) is equivalent to 

inr*(G) 



idiv(B^W^) = 




02 G 



B^W -N 

B^W -N = -- , ^ 

W{Q,Z2) = W{L,Z2) 
ly-iG) ^ I-^^gI dzidz2 = 



on B2(G) 
on Bi{G) 

on Bo(G) 



(A.4) 



where 



iB^W)iz,,Z2) 



-/"(Zl)22 i^(2l,Z2)+^ (l + (22F'(2l))2) §g{z,,Z2) ) ' 

That is, we have that X{G) is the solution of (A.l) in Y*{G) if and only if W{G) ~ X{G) o satisfies 
equation ([XI]) in r*(G). 



Moreover, if we define 



U{G) = W{G) 



\Y*{G)\ Jy>(g) 



W{G) 
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then, U{G) is the unique solution of 



f -ldWiB^U) = 



B^U ■ N 



in Y*{G) 
on B2{G) 
on Bi{G) 



^1 + (92G)' 

f/(0, Z2) = t/(i, ^2) onBo(G) 
/y*(G) U dzidz2 = 



(A.5) 



Also, we have that the variational formulation of problem (A.5) is given by the bilinear form 
P(5 : H{G) X H{G) ^ R 



{U,V) ^ 



y*(G) 



B^U -V 



dzidz2 



where H{G) is given by 

H{G) = {UeH\Y*{G))\U{0,z2) = U{L,z2)onBo{G), [ U = 0} 

Jy'(g) 



with the H^{Y*{G)) norm. Hence, we have that X{G) is the solution of (A.ll in Y*{G) if and only if 



mG)^XiG)oL,-^J^^^^XiG)oL, 



satisfies 



G' 



V 



dS, for all V G H{G). 



(A.6) 



(A.7) 



Observe also that the bilinear form associated to problem ( A.l I in Y*{G) is given by 

PG ■■ H{G) X H{G) ^ M 



and the weak formulation is given by 

Pg{U,V) = - 



{U,V) ^ 



G' 



VU ■ VV dzidz2 



Y*{G) 



VdS, ioi alW e H{G). 



,(G) ^T+wr 

As a matter of fact, we will be able to show the following 

Proposition A.l. Let us consider the family of admissible functions G G A{M) for some constant M , 
where A{M) is defined in \A. 



Then, for each e > there exists (5 > such that ifG,G£ A(M) with \\G — G\\ci{M) 5; then 

||f/(G)-C/(G)|U(y.(G)) <e. (A., 



Moreover, we obtain from (A. 8) 
where 



r{G) 



\r{G)~r{G)\<e 
dX{G) 



\Y*iG)\ 



(A.9) 
(A. 10) 



and similarly for r{G). 
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Proof. Obviously (A.9) follows straightforward from (A.8), the relation between U and X stated in (A.6), 



the definition of r{G) in (A. 10 1 and the fact that G and G are close in the C^-metric, and in particular in 
the uniform norm. 



Hence, we need to show (A.8 1. For this let pa and be bilinear forms associated to the variational 
formulation of (A. 5 1 for G and G respectively. We can get the following estimate 

\p^{U,V) - pg{U,V)\ 



< 



Iy'(g) ^ 



dzi \F 



< 



F dz2 \F 

F-1\\l 



dzi 
)dz2 



dzo' 



Go 
Gi 

\\FJl 



dV 




dU 


+ 


111 


dzi 


L2 


dzi 


L2 





dU 
dzi ) 
-^G\\\F' 



dz\dz2 



(^2 



dV 




dU 


dz2 


L2 


dz2 



dV 




dU 




dV 




dU 




dzi 


L2 


dz2 


+ 

L2 


dz2 


L2 


dzi 





ll-Fl 



dU 
dzi 



■G^\\F'\ 



dU 
dz2 



L2 



V 



L2 



(A.ll) 



Since 



F{zi) - 1 = and 

G{zi) 



G 



'(^1) (g(zi) - G(^i)) + G(zi) (g'(zi) - G'(zi) 
GWf 



we have 



11^ - 1||l<»((o,l)) < 7:^||G - G||loo((o,l)) and 
\F'\\l-°{[q.l)) < 773- ||G - G||ioo((o,L)) + Gi ||G' - G'||ioo(o,i,)j 



Then, due to ( A.ll[ ), we get 

\pq{U, V) - pg{U, V)\ < G (||G - G\\l^ + ||G' - G'||i^) \\U\\h\\V\\h, (A.12) 

where the constant G depends on M, Gq and Gi and therefore it can be chosen the same constant for all 
G,G G A{M). 

Now, since pc is a coercive bilinear form in H{G) (observe that we have imposed the condition /y,(-Q-) U = 

in F(G)) then, there will exist a constant c> such that Pg{U{G) - uIg),U{G) - U(G)) > c\\U{G) - 
t/(G)||^(.pj. But, to simplify, if we denote hy U — U{G) and U = U{G), we will have 

c||f^ - ^llff(G) < PGiU ~U,U-U)^ pg{U, U-U)- pg{U, U-U) 
< \pg{U, U^U)- p^{U, U-U)\ + \p^iU, U^U)- pg{U, u - U)\ 

G' 



< 



Bi{G) 



a 



1 

2 F 



\U -U\dS + C\\G~ G\\c^^)\\U\\h(g)\\U - UWUiG)- 



But using appropriate trace theorems in Bi{G), the fact that G, G £ A{M) so that they are uniformly 



bounded in G"'^(M) (and therefore 



G' 



G' 



only on M), then, we easily get that c\\U — U\\h{g) < G||G — G||c 



< G||G — G||c'i(R) for a constant G depending 
where we have used that we have 
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a uniform bound of for all G G A{M)^ which is easily obtained from the variational formulation 

(A. 7 1 and the fact that G E A{M). This shows (A. 8 1 and we conclude the proof of the result. □ 
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